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Adjoint Sensitivity Determination for 
Nonlinear Circuit Models 

Technical Field 

[0001 ] The present invention relates to computerized modeling of electronic 

circuits and, in particular, to computerized analysis of electronic circuit 
sensitivities by an adjoint method. 

Summary 

[0002] Electronic circuit computer simulators such as SPICE, HSPICE, Spectre, 

etc., are commonly used to model various characteristics of electronic circuit 
operation. These simulators formulate and solve the nonlinear algebraic 
differential equations associated with an electronic circuit design, as is known in 
the art. To improve electronic circuit design efficiency, the sensitivity of an 
electronic circuit to variations in circuit components is sometimes analyzed to 
identify components that are particularly sensitive to variations. In this regard, 
circuit sensitivity quantizes the effect on the performance of a circuit caused by 
some variation in a circuit component. 

[0003] Two main methods have been used for sensitivity computation: the direct 

method and the adjoint method. The direct method is efficient in the 
computation of several outputs with respect to one component. The adjoint 
method is efficient in the computation of the sensitivity of one output with 
respect to all the components. As a result, the adjoint method is typically more 
efficient in circuit sensitivity computations because most circuits have far more 
circuit components than outputs. A simple extension allows the adjoint method 
to generate the sensitivity of a function of several outputs. The adjoint method 
typically entails a first or original simulation of the circuit in its original form, 
determination of an adjoint circuit corresponding to the original circuit, and then 
a simulation of the adjoint circuit. 



1 



1051-006 

[0004] As is known in the art, nonlinearities in the original simulation may be 

handled by linearizing about a certain point. The quiescent or operating point is 
found by successively finer approximations about this point. Typically, after 
the operating point is found, the various nonlinear elements are treated as linear 
and are assigned the impedance values, etc that the elements would have at the 
operating point. This assures the linear circuit mimics the original nonlinear 
circuit at that point. 

[0005] The adjoint method is invoked in order to calculate the sensitivities of 

observables with respect to various resistors, etc. In the usual implementation, 
the nonlinear devices are fixed at their operating point values. But if the circuit 
is changed slightly, such as by adjusting the parameters, the nonlinear devices 
would in reality change their impedances, etc. as the circuit changes. This 
change is not reflected in the fixed values used in the typical adjoint method. 
This is sometimes referred to as the operating point shift problem. 

[0006] The present invention enhances the adjoint network method. The effects 

of nonlinear circuit elements are represented by augmenting the elements of the 
adjoint network. In particular, deviations away from linearity are represented in 
the original circuit by "fictitious" voltage sources. These voltage sources will 
map into "fictitious" current sources in the adjoint network. These sources are 
not static; they are directly proportional to the adjoint current through the 
branch corresponding to the nonlinear element. As such they may be classified 
as current-controlled current sources and are sometimes referred to as 
"correction" sources. 

[0007] Accordingly, an electronic circuit sensitivity analysis method for 

analyzing sensitivity of an electronic circuit model as represented by electronic 
circuit model data includes conducting a first computer simulation of the 
electronic circuit model and receiving results of the first simulation. A 
nonlinear circuit element is identified in the electronic circuit model and a 
nonlinear effect of the nonlinear circuit element is represented by applying a 
corresponding voltage source to the electronic circuit model. 
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[0008] An adjoint of the electronic circuit model is generated based upon the 

results of the first simulation, including mapping the corresponding voltage 
source into a current source in the adjoint. A simulation of the adjoint of the 
electronic circuit model is conducted and a circuit sensitivity analysis of the 
electronic circuit model is conducted based upon the results of the simulations 
of the electronic circuit model and the adjoint to it. 

[0009] Additional objects and advantages of the present invention will be 

apparent from the detailed description of the preferred embodiment thereof, 
which proceeds with reference to the accompanying drawings. 
Brief Description of the Drawings 

[001 0] Fig. 1 illustrates an operating environment for an embodiment of the 

present invention. 

[001 1 ] Fig. 2 is a functional block diagram illustrating functions performed by an 

electronic circuit optimization modeling software engine. 
[001 2] Fig. 3 is a component block diagram illustrating an implementation of a 

sensitivity analysis software for performing sensitivity analysis. 
[001 3] Fig. 4 is a circuit schematic diagram of circuit components of a simple 

operational amplifier circuit to illustrate typical electronic circuit modeling. 
[001 4] Fig. 5 shows an ideal voltage-controlled resistor to illustrate operation of 

an adjoint method. 

[001 5] Fig. 6 is a graph showing differences between exact and uncorrected 

adjoint circuit characterizations. 
[001 6] Fig. 7 is a graph illustrating an optimization error that can arise from 

uncorrected adjoint circuit characterizations. 
[001 7] Fig. 8 illustrates a nonlinear element in which the current across the 

element is a function of the voltages at the two terminal nodes of the element. 
[001 8] Fig. 9 illustrates the nonlinear element of Fig. 8 with corresponding 

substitute voltage sources for nonlinear effects and a corresponding adjoint 

circuit. 

[001 9] Fig. 10 is a graph of current-voltage relationships for a gain- specified 

voltage-controlled current source. 
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[0020] 
[0021] 



Fig. 11 is a flow diagram illustrating an adjoint sensitivity method. 
Fig. 12 is a schematic diagram of a circuit that is a generalization of the 



[0022] 



sample circuit of Fig. 5. 

Fig. 13 is a diagram of an adjoint circuit corresponding to the circuit of 



[0023] 



Fig. 12. 

Fig. 14 shows a multistage voltage-controlled resistor circuit that was used 



[0024] 



to test the adjoint network corrections described above. 

Fig. 15 is graph of input voltage versus operating point voltage for the 



multistage voltage-controlled resistor circuit of Fig. 14. 

[0025] Fig. 16 is a graph illustrating corrected and uncorrected sensitivities for 

multistage voltage-controlled resistors. 

Detailed Description of Preferred Embodiments 

[0026] Fig. 1 illustrates an operating environment for an embodiment of the 

present invention as a computer system 20 with a computer 22 that comprises at 
least one high speed processing unit (CPU) 24 in conjunction with a memory 
system 26, an input device 28, and an output device 30. These elements are 
interconnected by at least one bus structure 32. 

[0027] The illustrated CPU 24 is of familiar design and includes an ALU 34 for 

performing computations, a collection of registers 36 for temporary storage of 
data and instructions, and a control unit 38 for controlling operation of the 
system 20. The CPU 24 may be a processor having any of a variety of 
architectures including Alpha from Digital, MIPS from MIPS Technology, 
NEC, IDT, Siemens, and others, x86 from Intel and others, including Cyrix, 
AMD, and Nexgen, and the PowerPC from IBM and Motorola. 

[0028] The memory system 26 generally includes high-speed main memory 40 in 

the form of a medium such as random access memory (RAM) and read only 
memory (ROM) semiconductor devices, and secondary storage 42 in the form of 
long term storage mediums such as floppy disks, hard disks, tape, CD-ROM, 
flash memory, etc. and other devices that store data using electrical, magnetic, 
optical or other recording media. The main memory 40 also can include video 
display memory for displaying images through a display device. Those skilled 
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in the art will recognize that the memory 26 can comprise a variety of alternative 
components having a variety of storage capacities. 

[0029] The input and output devices 28 and 30 also are familiar. The input 

device 28 can comprise a keyboard, a mouse, a physical transducer (e.g., a 
microphone), etc. In addition, input device 28 includes an optical scanner that 
optically scans printed and other written documents or materials (together 
referred to as printed documents) to generate digitized images of them. The 
output device 30 can comprise a display, a printer, a transducer (e.g., a speaker), 
etc. Some devices, such as a network interface or a modem, can be used as 
input and/or output devices. 

[0030] As is familiar to those skilled in the art, the computer system 20 further 

includes an operating system and at least one application program. The 
operating system is the set of software which controls the computer system 
operation and the allocation of resources. The application program is the set of 
software that performs a task desired by the user, using computer resources 
made available through the operating system. Both are resident in the illustrated 
memory system 26. 

[0031 ] In accordance with the practices of persons skilled in the art of computer 

programming, the present invention is described below with reference to acts 
and symbolic representations of operations that are performed by computer 
system 20, unless indicated otherwise. Such acts and operations are sometimes 
referred to as being computer-executed and may be associated with the 
operating system or the application program as appropriate. It will be 
appreciated that the acts and symbolically represented operations include the 
manipulation by the CPU 24 of electrical signals representing data bits which 
causes a resulting transformation or reduction of the electrical signal 
representation, and the maintenance of data bits at memory locations in memory 
system 26 to thereby reconfigure or otherwise alter the computer system's 
operation, as well as other processing of signals. The memory locations where 
data bits are maintained are physical locations that have particular electrical, 
magnetic, or optical properties corresponding to the data bits. 
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[0032] Fig. 2 is a functional block diagram illustrating functions performed by an 

electronic circuit optimization modeling software engine 100 that is separate 
from an electronic circuit simulator 102. Simulator 102 may be any widely 
available simulator such as such as SPICE, HSPICE, Spectre, etc., or any other 
custom simulator. These simulators formulate and solve the nonlinear algebraic 
differential equations associated with an electronic circuit design. 

[0033] Optimization software 100 performs a sensitivity analysis 104 that 

generates sensitivity data that may be used for design optimization 106 and a 
mismatch analysis 108. Sensitivity analysis 104 provides information about 
how each circuit parameter affects the circuit output performance. This allows a 
circuit designer to identify parameter changes that will optimize output 
performance for key specifications. 

[0034] Design optimization 106 may be characterized as defining a a mapping of 
a multi-dimensional space in which each dimension corresponds to a different 
design specification into a one-dimensional space by means of an objective 
function of the circuit. Often the full dimensionality of the parameter space is 
not allowed, so the problem is further constrained by additional constraint 
functions. In this characterization, design optimization 106 amounts to adjusting 
circuit performance so that quality is optimized while necessary constraints are 
maintained. Design optimization 106 is completed when the parameters fall 
within the allowed constraints and quality is optimized. 

[0035] Mismatch analysis 108 minimizes the effect that individual components 

have on overall performance. Circuit components that match each other help to 
achieve the maximum circuit performance and reduce the risk that 
manufacturing process variations will cause the circuit to fail production tests or 
in an end-user system. Circuit components are considered to match each other 
when a circuit behavior is a function of the matching components' ratio, rather 
than the parameters of the individual components. 

[0036] The sensitivity data generated by sensitivity analysis 104 may also be used 

to conduct a statistical analysis 110 that is used in design centering 112. Design 
centering 112 is directed to optimizing manufacturing yields. Statistical analysis 
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110 applies process variation data 114 that represents circuit performance 
characteristics associated with particular manufacturing processes to the 
sensitivity data. In one implementation, statistical analysis 110 employs a 
Monte Carlo analysis that randomly varies every parameter in the circuit design, 
with "trials" being generated for each set of parameter values. In other 
implementations, statistical analysis 110 may employ Root Sum Square 
Analysis (RSS) or Worst Case Analysis (WCA, sometimes referred to as 
Extreme Value Analysis or EVA), as are known in the art. 

[0037] Fig. 3 is a component block diagram illustrating sensitivity analysis 

software 120 for performing sensitivity analysis 104. Sensitivity analysis 
software 120 includes a sensitivity software operation engine 122, sometimes 
called OpSens 122, which manages the operations for computing circuit output 
sensitivity to variations in circuit parameters and components. Sensitivity 
software engine 122 communicates with a separate circuit simulator 102 via a 
front-end interface 126 and a back-end interface 128. 

[0038] Sensitivity software engine 122 and circuit simulator 102 typically would 

employ different data formats for electronic circuit specifications and 
simulations. Front-end interface 126 and back-end interface 128 provide 
communication and data conversion between sensitivity software engine 122 
and circuit simulator 102. Simulator 102 provides circuit definition and 
simulation data to sensitivity engine 122 via back-end interface 128, and 
sensitivity engine 122 passes simulation commands to simulator 102 via front- 
end interface 126. Front-end interface 126 and back-end interface 128 are 
particularly adapted to the data format of circuit simulator 102 and allow 
sensitivity software engine 122 to be generically used with different circuit 
simulators having different data formats. 

[0039] For purposes of illustration, Fig. 4 is a circuit schematic diagram of circuit 

components of a simple operational amplifier circuit. In one implementation, 
the circuit definition and simulation data may be represented as a "netlist " 
which is a description of an integrated circuit design as is known in the art. In 
one implementation, a netlist description may have the following data structure: 
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*VDD0 21 5.0 
Vs0 21 5.0 
Vd20 0 5.0 

Vin in 0 AC 1 DC 1 PWL 0,-10 100u,10 

*Vin in 0 AC 1 DC 1 sin(0 3 lMeg) 

*Vin in 0 pulse(0 1 0 .In .In lus 2us) 

Rl inVN 10ktcl=-.2 

R2VNVO 10k tc 1=0.2 

VPVPODCO 

Xopl VO VP VN 20 21 amp 

.subckt amp VO VP VN NET67 NET34 

R77NET32VO 2E3 M=1.0 tcl=-.2 

C76 NET48 NET32 1E-12M=1.0 

*V70 0NET34 5.0 

*V12NET67 0 5.0 

M6 NET48 VP NET44 NET67 PMOS L=4E-6 W=30E-6 M=1.0 
M3 NET35 VN NET44 NET67 PMOS L=4E-6 W=30E-6 M=1.0 
M9 VO NET48 NET34 NET34 NMOS L=3E-6 W=154.2E-6 M=1.0 
M4 NET35 NET35 NET34 NET34 NMOS L=4E-6 W=15E-6 M=1.0 
M7 NET48 NET35 NET34 NET34 NMOS L=4E-6 W=15E-6 M=1.0 
M2 NET54 NET54 NET34 NET34 NMOS L=32E-6 W=3E-6 M=1.0 
M8 VO NET54 NET67 NET67 PMOS L=4E-6 W=200E-6 M=1.0 
Ml NET54 NET54 NET67 NET67 PMOS L=4E-6 W=12E-6 M=1.0 
M5 NET44 NET54 NET67 NET67 PMOS L=4E-6 W=30E-6 M=1.0 
.ends amp 
.end 



[0040] The circuit definition and simulation data are stored in a simulation results 

database 130. Sensitivity engine 122 cooperates with a sensitivity calculator 
132 to determine sensitivity data from the circuit definition and simulation data 
in simulation results database 130. Sensitivity calculator 132 may be controlled 
by calculator scripts 134. The sensitivity data are held in a sensitivity file 136 
that may be used, for example, in statistical analyses 110. 

[0041] Sensitivity engine 122 uses circuit simulator 102 to simulate the selected 
circuit and, using data gathered from this nominal simulation, builds one or 
more modified adjoint circuits. Sensitivity engine 122 then uses circuit 
simulator 102 to simulate the one or more adjoint circuits. The outputs from 
these nominal and adjoint simulations are used to populate database 130 with 
sensitivity information for each variable in the circuit network, which allows 
identification of the circuit components to which the circuit outputs are most 
sensitive. This sensitivity information can be used for design optimization, 
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centering, and other analyses. For example, sorting routines can assist in 
identifying the most sensitive parameters, thus simplifying the task for designers 
to reduce circuit sensitivity to component variations. 
[0042] The sensitivity is a measure of the effect on circuit performance P 

(sometimes called an observable) due to the variation of some circuit element x: 
x dP 

S p = . Sensitivities relating to multiple parameters, as is typical, may be 

x P dx 

x dP 

expressed as S% = — where i indexes the parameter. The sensitivity is an 

P dx t 

indispensable tool for the design, test, and optimization of circuits. For 
instance, a large class of optimization methods rely heavily on gradients, which 
dP 

are denoted by — and give information essentially equivalent to sensitivities. 
dx 

An efficient and accurate method of generating sensitivities for circuits is thus a 
crucial need for circuit designers. 

[0043] The qualifier "efficient" is important. Sensitivities may be generated by 
numerical differentiation, (sometimes called the brute force, or the direct 
method). With this method, one parameter at a time out of the N is varied 
slightly, or perturbed, the simulation is rerun , and P is calculated again. The 
numerical derivative, and hence the sensitivity, is constructed from the 
difference between this value of P and the original. This is quite 
straightforward, but with N adjustable parameters it requires N+l simulations at 
every step. With simulations correspondingly slow (likely on the order of order 
0(N) time at least), this results in at least 0(N 2 ) complexity. With elements 
running into the thousands, unmanageable run times are quickly reached. 

[0044] The adjoint method, or the adjoint network method, uses techniques of 
linear algebra to map a circuit into a related "adjoint" circuit. A single 
simulation of this adjoint circuit yields the derivatives of any circuit 
performance measure (observable) P with respect to all parameters of interest: 

dP 

^,/ = (l,2,..JV). 
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[0045] This is a very efficient method for two reasons: (1) After the original 
simulation, only one simulation per performance measure ( P ) is needed, 
regardless of the number of adjustable parameters. (2) The adjoint circuit is 
linear, even if the original circuit is not, so adjoint simulations are likely to be 
faster than the original. However, the statement (3) above points out a problem 
with the conventional adjoint circuit method. It can sometimes give inaccurate 
and misleading results when nonlinear elements are present in the original 
circuit. 

Analytical Example: Conventional Adjoint Method with Nonlinear Elements 

[0046] Fig. 5 illustrates an ideal voltage-controlled resistor. The resistance R s is 
fixed at a factor II a times the voltage at the first node of the resistor v x . This 
problem may be treated analytically. There is one circuit equation 

(1) v-I(R 1 +R 2 +R s )^0 
and one equation for the resistor 

(2) v-IR x =aR s 
so the solution is 

-(aR x +aR 2 - v) + ^(a^ + aR 2 - v) 2 + 4aR 0 v 

(3) R s = - 

from which the other quantities may be easily determined. 

[0047] The conventional implementation of the adjoint network method 

substitutes the above value of R s back into the circuit equation. (In the more 
general nonanalytical case, this value would have been determined by a 
simulation.) This value of R s is fixed. Given the values (v,R l9 R 2 ) , the JR 4 
value chosen above will be consistent, and the circuit might as well be linear. 

[0048] Consider what happens when one of the parameters, say R x , changes 

slightly, with the voltage v 2 chosen as the observable. By an elementary 
application of the adjoint network method (or just by direct differentiation) the 
following may be derived: 
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(4) 



dv 2 
dR x 



-vR 2 _ -vR 2 



(R x +R 2 +R s ) 2 r2 



where R = R x + R 2 + R s . The U subscript designates "Uncorrected", as explained 
below. By differentiating equations (1) and (2) above with respect to R x , and 
dR 

eliminating — - , the exact, analytical result is: 
dR, 



(5) 



2 



dR x 



dl n la -I 



2 



= R 2 — = R 
2 dR x 



IR X ~aR 



where the E subscript means "Exact." The Exact and Uncorrected quantities are 
not the same, as the plot of Fig. 6 shows. This discrepancy can affect 
optimization. 

[0049] Consider an optimization that holds R 2 , a, and v constant, and maximizes 
the power I 2 R 1 across R x . Fig. 7 is a graph of I 2 R X vs. R x with 
R2 = 1000, v = l,a = .01 . This graph shows that that the maximum power 
occurs at about /? 1 =900. But if the uncorrected values had been used (in other 



dR s _ „._™„ 1JU ^ n_ r2 D _ v 2 i? x 
dR x 



words set z^r = 0), we would have P = I 2 R X = — ^ and 



R 

= 2Rl ^ with zero at R~+R=R,. This relation, as well as the 

dR x - 2 - 3 
1 R R 

equation v-IR x =aR s , is satisfied at the point R 2 = 1000, R s = 500, R x = 1500 , 
/ = 1/(3000) . However, this point is clearly not the true power maximum. 
(Notice this is not a problem when numerical derivatives are computed. There, 
slight variations of elements will be countered by self-adjustment of the 
nonlinear devices, and when the simulator converges the observed values are 
the correct ones.) 
Operating Point Sensitivity Correction due to Nonlinearity 
[0050] The present invention enhances the adjoint network method. The effects 

of nonlinear circuit elements are represented by augmenting the elements of the 
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adjoint network. In particular, deviations away from linearity are represented in 
the original circuit by "fictitious" voltage sources. These voltage sources will 
map into "fictitious" current sources in the adjoint network. These sources are 
not static; they are directly proportional to the adjoint current through the 
branch corresponding to the nonlinear element. As such they may be classified 
as current-controlled current sources and are sometimes referred to as 
"correction" sources. 

[0051 ] As might be expected, the coefficient by which the correction source is 

proportional to the adjoint current is itself a multiple of the nonlinearity of the 
device, as measured by the derivative of the impedance or admittance with 
respect to the voltage across the device. Thus, linear devices require no 
correction. Several common descriptions of nonlinear elements are illustrated 
with respect to: 

(1) admittance 

(2) impedance 

(3) gain ( g ) characterization of voltage-controlled current source 
Although the solution to the operating or quiescent point problem is 
independent of the description, it is instructive to tailor the algorithm to suit the 
representation at hand. 

Admittance case 

[0052] Fig. 8 illustrates a nonlinear element in which the current across the 

element is a function of the voltages at the two nodes. For most cases of interest 
the current is a function only of the difference between the voltages, v 2 -v l . In 
addition, one of the nodes will be at ground. However, these conditions do not 
always hold; even in the simple voltage-controlled resistance above, for 
instance, they are violated. So the most general case will be solved, where 

/ = i(y v ) . It is useful to define an admittance y = — - — where it is 

v 2 -v x 

understood that y = y(v 15 v 2 ) , unlike the linear case. This definition will allow 
the effect of the nonlinearity to be isolated. The case of one control will first be 
described, and the case of multiple controls will be described later. 
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[0053] Fig. 9 illustrates the nonlinear element of Fig. 8 with its state at the 

quiescent or operating point O given by the values (v 10 , ,I Q7 y 0 = — — — ) . 
For perturbations about this operating point O, the admittance will be given by: 

y = yo + "^" Av i +^-Av 9 , where &v x =Vj -v 10 ,Av 2 =v 2 -v 20 . 
As illustrated in Fig. 9, the nonlinear element may be "replaced," to first order 
in Av 19 Av 2 , with a linear element of nominal admittance y 0 , plus two voltage 

sources (J l kv l and <J 2 Av 2 where: 



(6) a x ^ 



2 



y 0 y 0 
<h> , \ ^y 



(7) <T 2 = — 

y 0 y 0 

[0054] The nonlinear element is replaced in that the currents of the nonlinear and 

linear elements are the same so that the branches are indistinguishable. The 
current in the nonlinear element before replacement: 

(y 0 +|^A Vl +^-Av 2 )(v 20 + Av 2 -v 10 -Av,) 2 ) 

can be demonstrated to be is equal to that after replacement: 

y 0 O20 + Av 2 - v io - Av i + Av i + a i Av > 
to first order. First, the zero-order terms are clearly equal; the first-order terms 
are: 

(|^Av x + -|^-Av 2 )(v 20 -v 10 ) + y 0 (Av 2 - AvJ 

and 

y 0 (Av 2 - Av x + cr,Av 1 + <7 2 Av 2 ) , 
respectively. Subtracting the v 0 (Av 2 - Av x ) from both, and using equations (6) 
and (7) these are trivially equal. 
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[0055] This shows that a nonlinear effect may be modeled by linear elements that 

can be easily handled by the adjoint process. To be precise, the nonlinear effect 
has been incorporated into a pair of voltage sources, cr l Av l = cr l (v 1 - v 10 ) and 
a 2 Av 2 = <T 2 (v 2 - v 20 ) . In going to the adjoint network, the constant terms in 

the sources disappear, as known in the art, and the voltage sources become 
corresponding current sources 

dy dy 

aj = ?\l 0 l and a 2 J=%/ 0 J 

Jo ?o 

A 

where the / is the current through the adjoint element corresponding to the 
nonlinear element as shown in Fig. 9. The desired derivatives will be products 

A 

like I 0 1 , as described below in greater detail with reference to a full adjoint 

network construction. 
Impedance case 

[0056] The analysis of the shifting operating point problem using admittances is 

solved above, but admittances may not necessarily be the most accessible 
quantities. If impedances are specified, the extension is straightforward. Since 

R = —, then: 

y 

3(— ) 

dv x y Q R dv x y 0 

(9) a x = -L and likewise <T 2 = -I 0 . 

[0057] Fig. 11 is a flow diagram illustrating an adjoint sensitivity method 150. 

Sensitivity method 150 includes a step 152 of conducting a first or original 

simulation of a nominal electronic circuit specification or model and receiving 

results of the first simulation. 
[0058] A step 154 represents the effects of nonlinear circuit elements (i.e., 

deviations away from linearity) in the original or nominal circuit by applying 

fictitious voltage sources. 
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[0059] A step 156 generates an adjoint of the nominal electronic circuit based 

upon the results of the first simulation, including mapping the fictitious voltage 
sources representing nonlinear circuit elements into fictitious current sources in 
the adjoint network. 

[0060] A step 158 conduct a simulation of the adjoint of the circuit by the 

electronic circuit simulator and receives results of the adjoint simulation. 

[0061] A step 160 calculates a circuit sensitivity analysis of the nominal 

electronic circuit specification based upon the results of the simulations of the 
nominal electronic circuit and its adjoint. 

Adjoint Network Derivation 

[0062] It is instructive to explore in more detail the changes introduced into the 

adjoint circuit. Although one could simply write down the adjoint terms 
corresponding to the added voltage sources, it is more instructive to re-derive 
the results from first principles. This derivation begins with a circuit (Fig. 12), 
which is a generalization of the sample circuit of Fig. 5, and uses impedances 
rather than admittances because of the topology of the circuit. R* designates a 
variable impedance of the nonlinear element. The circuit equations are 
(11a) v-IR^v, =0 

(lib) v.-IR'-v^O 
(11c) v 2 -IR 2 =0 
and the nonlinear device equation is 

(lid) R*-f(y x ,v 2 ) = 0 
[0063] Consider measurements of v 2 , meaning that the following are desired 

dv 2 dv 2 fo 2 
dR^ ' dR 2 ' dv 

The Lagrangian multiplier formulation is known in the art. R. A. Rohrer, "Fully 
Automated Network Design by Digital Computer: Preliminary Considerations," 
Proa IEEE, 55, 1929-1939 (1967). We construct the form consisting of the 
system variable v 2 whose derivatives are of interest, combined with the circuit 
constraint equations (11): 
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(12) 

r = v 2 + 4 (v - IR, - v t ) + /l 2 (v : - IR* - v 2 ) + A 3 (v 2 - ZR 2 ) + /l 4 (i? * - / (v, , v 2 )) 
Then the equations are solved for the system variables: 



[ dv 1 



— = ~\ -i 4 



= o 



Sv : 

sr 



. — 1 — X 2 + yl 3 — A 4 



-0 



Si 

where the | designates that values are to be taken at the operating point. 
[0064] Ordinarily this would complete the analysis, however, R* is now a system 

variable: 

SR 

The factor /t 4 is eliminated by the last equation so that the adjoint circuit (Fig. 13) 
may be read off from the first three equations. These equations require the 

introduction of two current sources, A 2 I -f- and A 2 I , in exact agreement 



with the earlier results. The following can also be expressed: 

dv 2 _ ST 
dR, ~ SR, 

ST 



(13a) ^r = -^- = -H 

• = -a, 



(13b) ^ 



(130 



dR 2 SR : 
ST 



= 4 



3v Sv 

Notice, as mentioned before, that one can retrieve the derivatives of v 2 with 
respect to all parameters of interest with just one adjoint simulation. 
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Comparison of analytical and adjoint network computations 
[0065] These results can be tested directly with the voltage-controlled resistor 

case in which f(y,,v 2 )=—,^- = —and-^- = 0. The A equations become 

a dv x a dv 2 

-A.R.-A^ -A 3 R 2 = 0 

/L = A+—A 2 =>A l = (l-—)A 2 
a a 

A, = 1 + A, => A, — A 2 — 1 



dv 

We wish to show — - = -I A so we solve for A 1 to get 
dR i 

-R 2 I(l--) R-l(a-I) 

(14) -AJ = ^ = 2 _ 

RJl-^ + R* +R 2 IR,-aR 
a 

which precisely matches the earlier result. 
Experimental Results 
[0066] Fig. 14 shows a multistage voltage-controlled resistor circuit that was used 
to test the adjoint network corrections described above. A plot of the voltages at 
various stages is shown in Fig. 15. 
[0067] Fig. 16 is a plot directed to V4, the last-stage voltage, and its sensitivity 

with respect to v, — (note the scale change). Here a geometric interpretation 
dv 

is possible. — should match the actual slope of the V4 curve . The corrected 
dv 

and uncorrected sensitivities are illustrated, and the disparity between the 
corrected and the uncorrected is clear. The accuracy of the corrected sensitivity 
(i.e., how well it matches the derivative) is also clear. As a check, numerical 
differentiation was also obtained by varying v slightly at each point. The 
numerical derivatives match precisely the corrected adjoint derivatives. 
[0068] The CPU requirements for this case were as follows: 
Original simulations: 1.68 seconds 
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Uncorrected Adjoint Network simulations: .87 seconds 
Corrected Adjoint Network simulations: .78 seconds 
Numerical derivatives (to check Corrected Adjoint): 2.85 seconds 
Simulations were performed on an in-house simulator based on relaxation 
methods. Surprisingly, the corrected adjoint method took slightly less time than 
the uncorrected method. This is probably not significant. Both the original 
simulation and the numerical differentiation took substantially longer. This 
example shows convincing evidence for the efficacy of the corrected adjoint 
network sensitivity method for nonlinear elements. 
Summary 

] The following table summarizes the adjoint correction mechanisms 
described above. This is specialized to the case where one of the nodes of the 
nonlinear device is at ground, the other at voltage v. Here 7 0 is the nominal, 

A 

operating point current through the device, 7 is the adjoint current through the 

A 

device, and 7+ is the additional adjoint current source correction. 
Table I: Characterization and Adjoint Circuit Correction for Nonlinear Devices 



Characterization of 
nonlinear device: 


I-V Relations satisfied: 


Adjoint Current Source 
Correction 


admittance y 


I 

y = - 

V 


dv y 0 


impedance R 


R = V - 
I 


/ -^7 7 



[0070] Automated design of analog and mixed-signal circuits is becoming more 
and more critical with increases in circuit size. Efficient and accurate 
calculation of sensitivities is a crucial part of that automation process. The 
adjoint network method is an efficient method for calculating sensitivities, but 
in its conventional form, is partially flawed when used in the presence of 
nonlinear circuit elements As described above, the method may be corrected 
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by augmenting the adjoint circuit with current-controlled current sources whose 
coefficients are proportional to the measure of the nonlinearity of the device. 

[0071 ] The operating point correction introduces current sources into the adjoint 
circuit. Each such source is a current-controlled current source, which is a 
standard circuit construct; however, an element of choice exists with regard to 
its implementation. Most simulators function by repeatedly solving the linear 
system which we indicate symbolically by = b where the A represents the 
passive elements of the circuit, b is a source vector, and x is the state vector. 
Even nonlinear circuits can be solved this way, by iteration; this is how the 
operating point of a nonlinear circuit is found. The point is, the current source 
we have constructed above is now a linear function of the state vector. We can 
represent this by b = Jx + v , say, where J is constant. Then we have a choice 
whether to solve the system Ax = Jx + b\ ov Ax- Jx = V The first can be 
solved iteratively; the second in one step if, as is typically the case with the 
adjoint circuit, A is constant. Adjoint simulations are typically much faster 
because they are linear; it may or may not be worthwhile to modify A in order 
to get a quicker solution. 

[0072] Having described and illustrated the principles of our invention with 

reference to an illustrated embodiment, it will be recognized that the illustrated 
embodiment can be modified in arrangement and detail without departing from 
such principles. In view of the many possible embodiments to which the 
principles of our invention may be applied, it should be recognized that the 
detailed embodiments are illustrative only and should not be taken as limiting 
the scope of our invention. Rather, I claim as my invention all such 
embodiments as may come within the scope and spirit of the following claims 
and equivalents thereto. 
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